Quantum geometric phase in Majorana's stellar representation: 
Mapping onto a many-body Aharonov-Bohm phase 
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The (Berry-Aharonov-Anandan) geometric phase acquired during a cyclic quantum evolution 
of finite-dimensional quantum systems is studied. It is shown that a pure quantum state in a 
(2,7 + l)-dimensional Hilbert space (or, equivalently, of a spin- J system) can be mapped onto the 
partition function of a gas of independent Dirac strings moving on a sphere and subject to the 
Coulomb repulsion of 2 J hxed test charges (the Majorana stars) characterizing the quantum state. 
The geometric phase may be viewed as the Aharonov-Bohm phase acquired by the Majorana stars 
as they move through the gas of Dirac strings. Expressions for the geometric connection and 
curvature, for the metric tensor, as well as for the multipole moments (dipole, quadrupole, etc.), are 
given in terms of the Majorana stars. Finally, the geometric formulation of the quantum dynamics 
is presented and its application to systems with exotic ordering such as spin nematics is outlined. 
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The concept of geometric phase associated with a 
cyclic quantum evolution [l], Q has by now become a 
central unifying concept of quantum mechanics • Its 
importance stems from the fact that its local expression, 
the geometric curvature, controls the quantum dynamics. 

The prototype of geometric phase is that of a spin 1/2: 
in this case, it can be thought of as the Aharonov-Bohm 
(AB) phase of a unit electrical charge moving on the 
Bloch sphere in the field of a Dirac monopole of unit mag- 
netic charge located at the center of the Bloch sphere [l| . 
For a spin- J system evolving within the manifold of co- 
herent states (CSs) (i.e., quasi-classical states), this inter- 
pretation goes over straightforwardly, the electric charge 
being now of magnitude 2 J For a spin- J system in an 
arbitrary quantum state, except for for the specific case 
of a quantum evolution consisting only of global rota- 
tions (he., due to a, possibly time-dependent, magnetic 
field) the geometric phase associated with an arbi- 
trary cyclic evolution of a spin J system Q does not 
seem to be amenable to any physically appealing AB-likc 
interpretation, which is somehow unsatisfactory. 

In the present Letter, I propose a novel theory of the 
geometric phase of spin systems, based upon a mapping 
onto a (fictitious) many-body system. From the latter 
emerges quite naturally a novel AB-like, physically trans- 
parent, understanding of the geometric phase of spin sys- 
tems. Since any system with a finite dimensional Hilbert 
space of dimension 2 J + 1 can be thought of as a spin- J 
system, the present study actually holds, at least for- 
mally, for any finite quantum system. 

Since any two state vectors \ i $>\ ) and \^2^) of the 
spin-J Hilbert space H"' = C 2J+1 — {0} satisfying 
1*60) = cl*^) (with ceC,c/0) yield the same ex- 
pectation value for any observable, they represent the 
same physical state and belong to the same equivalence 



class (|^2 )~l^ r i''')); thus the manifold of physical states 



,60 



(projective Hilbert space) is the quotient space of equiv- 
alence classes 7>60 = U {J) /~ = CP 2J [2\. As in Ref. 0, 
the approach to be used here is purely geometric and re- 
lies upon Majorana's stellar representation Q for T>( J \ 
Majorana's representation is most easily understood by 
noticing that spin-J states can be obtained as fully sym- 
metrized states of system of 2 J spins 1/2 0. This idea is 
at the heart of the Schwinger boson (SB) representation 
in which the spin-1/2 CS pointing along the direc- 
tion n of spherical angles and ip is In' 1 / 2 ') = at |0), 
with at = cos (f ) at + sin (|) e 1¥, a|. Let us pick 2 J (non 
necessarily distinct) unit vectors {ui, . . . , U2,/} 

form the state |W 



(J)\ - 
u / — 



(llf'.<< n 



U, and 
. Obvi- 



ously, being a superposition of states with 2 J SBs, such 
a state is a spin-J state. In particular, the states |fi60) = 



r(4: 



2J| 



uct is given by 



are the spin-J CSs their scalar prod- 
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,iJE(z,n 1 ,n 2 ) 



where S(z,n 1 ,n 2 ) is the oriented area of the spheri- 
cal triangle (z, 111,112), and they satisfy the following 
resolution of unity: lj = 2^±1 J d 2 n |n60)(n60|. 
The rotated SB creation and annihilation operators sat- 
isfy the commutation relations: = [afi,ajy] = 
and [a fi ,at,] = (n (1 / 2) \h"^/ 2) ) . Let us introduce 
the CS representation *u = (n (J) 1 ) , which 
is a wavefunction over the sphere S 2 , with probabil- 
ity distribution Q^\h) = \^^\h)\ 2 (Husimi func- 
tion). Simple algebraic manipulations yield ^/^(h) = 

Iliii (Hr^ e is ^ ft <- a ')) 1/2 . Conversely, using the de- 
composition in the familiar \JM) basis, one sees that 
a generic spin-J state vector |\E'6')) can \, e expressed as 
I 'J 60) = P^(at,at)|0), where P*(a|,aj) is a homoge- 
nous polynomial of degree 2 J of at and at, which can 
be factorized (up to an unimportant prefactor) in the 
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2 J-body classical system 14]. The partition function can 
be expressed in terms of the pairwise chordal distances 
between the MSs, 2d i? = 1 — iij • u } -, as [l5l 
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FIG. 1. Typical diagrams used to compute fl{J (a), D 



,(J,n) 

(b) , and i?u™i/ ( c ) ■ The solid dots represent the MSs (labeled 
from 1 to 2J), the open dots (labeled fi, u, ...) represent 
the auxiliary stars used to compute the multipole moments. 
The rules are: (i) draw all possible distinct diagrams with n 
pairing links (each dot, solid or open, may be linked only once 
in a given diagram); (ii) calculate the contribution of each 
diagram as indicated below, and then sum over all diagrams; 
(iii) an unlinked solid dot yields a factor 1; (iv) an unlinked 
open dot yields a factor 0; (v) a link between 2 solid dots i 
and j yields a factor dij ; (vi) a link between a solid dot i and 
an open dot /i yields a factor iii M ; (vii) a link between 2 open 
dots p, and v yields a factor — 2<5 M „. The diagrams (a), (b), 

(c) shown here yield the contributions diidis, £5^23^46, and 
U3,j.U(iv to -Dy' 2 ', Dfi'® and Dfi'*£, respectively. 



above form ) w ith a unique multiset U (Fundamen- 
tal Theorem of Algebra). Thus, V {J) can be univocally 
parameterized by the constellation U of 2J Majorana 
stars (MSs) (zeros iii of the Husimi function), and 
can be taken as a fiducial state in T-L^ to describe V™. 
Addition or removal of n stars to/from a given spin- 
J constellation generates a state of spin (J + n/2) or 
(J — n/2), respectively. While the SB formalism was in- 
spired by Majorana's representation (l2j |. the underlying 
geometric aspects have not been fully explored so far; to 
carry out this task is one of the aims of the present paper. 

Noticing that the two-dimensional (2D) Coulomb in- 
teraction on the sphere is V"(ui,u 2 ) = — ln(di 2 ), where 
2rfi2 = 2 sin 2 (6*12/2) = 1 — Ui • U2 is the chordal dis- 
tance between Ui and U2 [IH, and writing Q^\ii) = 
ex P (-/3l\j(n)), with = 1 and U v (h) = Y.Zi ^( fi , Ui), 
we can interpret the (rescaled) norm of l^u'), 



with [J] = J (resp. [J] = J — 5) for 23 even (resp. 
2 J odd), and where the expression of £)£'"' in terms of 
diagrams is explained in Fig. [TJ 

Let us now show how the expectation value of the var- 
ious multipole moments can be expressed in terms of the 
MS. They are obtained from the expectation values of 
the irreducible spherical tensor operators y, m (J), which 
in turn have the following P-representation 16 1 



3T(J)= 



(2J + 1 + Q1 2J + 1 
(2J+ 1)!2 ; I^T" 



d 2 n|n^)Y7"(n)(n 



W I 



(3) 

For the expectation values of the dipole moment and 
the quadrupole moment = J ^ J "+ J " J " _ 
(where fx, v, . . . label the cartesian axes), this yields 



(J fl ) = (J + l)(n fl ), 



where 



/ s2 d 2 fi/(n) Q#(fi) 



J s2 d 2 nQg ) (n 



(4a) 

S -f) , (4b) 
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and (ri^n,, 



I remark that 



To calculate the averages (n p 
if we form the spin-( J + 1/2) state U' obtained from the 
spin- J state U by "adding" the star u', and the spin- 
( J + 1) state U" obtained by adding one further star u", 
the corresponding fictitious free energies are given by 



F(U') 
F(U") 



F(U)+F(l/2)-ln[l-^(rV)] (6a) 
F(U) + 2F(l/2) 

In [l - (vf + u'fiinp) + u^u^h^hv)] ,(6b) 



Z(U) = 



2J + 1 



1 



1 

47T 



exp 



(1) 



as the (fictitious) partition function, at inverse tempera- 
ture P = 1, of a classical gas of independent particles (of 
density Q^(n)) living on the sphere and interacting via 
the 2D spherical Coulomb repulsion with 2 J fixed test 
charges located at the Uj's. The corresponding fictitious 
free energy, F(U) = — /3 _1 lnZ(U), expresses a fictitious 
indirect interaction among the MSs, mediated by the gas 
particles in thermal equilibrium at temperature /3 _1 ; we 
thus obtain a mapping of a spin- J quantum state onto a 



where F(l/2) = In 2 (here and further below, Einstein's 
convention of summation over repeated is used, unless 
explicitly specified) . Thus we see that we can obtain the 
dipole and quadrupole moments by adding 1 or 2 auxil- 
iary stars, respectively, from the variation of the free en- 
ergy as these auxiliary stars are moved around the sphere. 
A careful but straightforward calculation yields: 



1 ES 1/21 (-l) n (2^+l-n)!^ n) 



2 ( J+1 ) El;io(-i)"(2J- 



EL J + 1I (-l)"(2J+2- ?v ^ u ^ 
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(7a) 



2(J+l)(2J+3) £M (-l)«(2J-n)!£>g 



(J,n) 



(7b) 
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where the expressions of D^'^ and D^'™^ in terms of 
diagrams are given in Fig. [TJ For example, the spin-1 
dipole and quadrupole moments are, respectively, 



(44 : 



2-d 12 ' 

1 (ux^UZv + U^Ulv 



2-d 



12 



(8a) 



-n r u 2 U -f). (8b) 



The extension of this procedure to higher-order multi- 
pole moments is straightforward. In turn, the method 
presented here allows to express, in terms of the MSs, 
the expectation value H(XJ) of the Hamiltonian and its 
derivatives 8 , which will be used further below to 
describe the quantum dynamics. 

Let us now come to the geometric phase and the quan- 
tum metric. The geometric phase acquired as the systems 
is parallel-transported along a closed circuit C in is 
given by [l|, [2| 



2.1 



A ■ dU = (p a, • du 



i=l 



(9a) 



2 J 2 

\ E E / f?fdu?AduP, (9b) 



with 



= 

I ij 

a/3 



da.? 



d^a? = -2 hn(htf), 



(10a) 
(10b) 
(10c) 



where 



ra 



(■/) 
u 



d (resp. a) indicates derivative of the b 

(resp. ket onr y- I n tne above equations, 

1, 2 label some spherical coordinates for the MS, 



(* 
a, 

with the tangent unit vectors ej and ef — x ej. In 
going from Eq. (|9aj) to Eq. (|9b[) . Stokes' theorem has 
been used, and S is an oriented surface bounded by the 
oriented path C. Here, af and are, respectively, 
the (gauge-dependent) Berry connection and the (gauge- 
independent) Berry curvature tensor; they have the phys- 
ical meaning of a "vector potential" and of a "flux den- 
sity", respectively, in V^ J \ The other important geomet- 
ric structure is the quantum metric (Fubini-Study met- 
ric), corresponding to a distance between |^) and |$) de- 
fined as D F s(*,*) = 2 arccos ( j^^^^m ) , whose 
infinitesimal expression is ds 2 = gf^dufdiij, with metric 

tensor gff = 4 Re(/i°f ) 0. 

The direct calculation of the geometric phase Q 
is complicated because of the need of taking care of 
the commutation relations among SB operators, yield- 
ing physically obscure results. This difficulty can be 
overcome by inserting the CS resolution of unity be- 
tween bra and kets, and, after some algebraic manip- 
ulations, one obtains a.; = (aj(n)), where a^(n) = 



"X f i-^"a- ~ l - fl"fl ) ^ s rea dily seen to be the vector po- 
tential, at iii, due to a (unit flux) Dirac string enter- 
ing the sphere along the z-axis and exiting at n. This 
means that the particles of our fictitious classical gas ac- 
tually carry a Dirac string; thus the MSs are surrounded 
by a flux density (of total flux equal to that of a Dirac 
monopole of unit magnetic charge) proportional to the 
gas density (n). The geometrical phase is then nat- 
urally interpreted as the AB phase acquired by the MSs 
as they perform a cyclic motion on the sphere. The 
most salient feature of this novel interpretation is the 
"fluid" character of the flux density, which results from 
the Coulomb repulsion between the flux carrying gas par- 
ticles and the MSs. For a CS circuit, Berry's result [1| is 
recovered, albeit with a different AB-like interpretation. 
Skipping technical algebraic details, the final expression 
for the Berry connection is (no Einstein convention here) 



-1 



Z X U, 



2 \ 1 - Z • U,; 
— 1 / ZXUj 

~2~ \ 1 - z • Ui 
-1 ( zxu, 
2 \ 1 - z • u, 



n 



1 n • Uj 

(n)' x u, ; 



x u. 



9u, F x u. 



(11a) 
(lib) 
(11c) 



In Eq. (jllbp . the notation (/(n))'j indicates that the av- 
erage is taken for the spin-(J— 1/2) state obtained by 
removing the star from the Majorana constellation of 
the spin- J state U; similarly (/(fi))y, to be used further 
below, indicates the average taken over the spin- (J— 1) 
obtained by removing the two stars and Uj, The first 
term in the above equations corresponds to the uniform 
flux density of a Dirac monopole, while the second one 
corresponds to the non-uniform part (with zero average) 
of the flux density. 

The metric and Berry curvature and tensors are ob- 
tained in a similar manner. For the former, one gets (no 
Einstein convention in Eqs. (ll2a|[T4"10 



a/3 - - 9 2 F 



9if = <%<W +e? 



2 F 



dxii ® duj 



= e?-d ij -e P +Je°'-d ij -Je?, 
with Jej = ef, Je 2 = -ej, and 



(12a) 
(12b) 



dij = 



n 



Ui • n 
n 



n 



u, n 



(12c) 
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For i j, this yields 
d. 



1 - (iii + Uj) • (&)'.. + u, • (n <g) n}^ • iij 



n 



(ft);- 



l-(n)' r u 4 l-(n);--u/ 



whereas for i = j, one gets 



1 - (WO 2 

Ci — <a>j - a*) 2 



(13) 



(14) 



The Berry curvature tensor is obtained from the metric 
tensor from the following identities 



2 ,12 - of 21 - a 11 - a 22 



»*3 



(15a) 
(15b) 



Mathematically, this follows from the Kahlerian nature 
of the projective Hilbert space [f| . 

To describe the quantum dynamics, I write down 
Schrodinger's equation for the state vector \tp(t)) = 

r , which yields tp = --ff(U) + A • U 



>> 



(*U(t)l*U(t)> 

(setting ft = 1); the latter result is nothing but as the 
infinitesimal version of the Aharonov-Anandan decom- 
position of the total phase variation into the dynami- 
cal and geometric terms 0. It is not gauge invariant 
with respect to a change of the phase choice for the 
fiducial states |\&u). In order to obtain a gauge in- 
variant equation of motion, we use the fact that, due 
to the unitarity of Hamiltonian evolution, the relative 
phase of any two (non-orthogonal) states, defined as 
<pi2(t) = axg{tpi(t) \tp2(t)), is time-independent. Doing 
this for states at U and U + 5XJ, one finally obtains 



(16) 



The above equation has the form of the classical equation 
of motion of a system of 2 J coupled particles evolving 
on a spherical phase space, with symplectic form given 
by hf^du? A diij, and a Hamilton function given by 
H(JJ) [18]. An alternative (equivalent) formulation of the 
spin dynamics in terms of the Majorana stars has been 
given earlier by Leboeuf [lil ]: however, the symplectic- 
Hamiltonian nature of the dynamics is displayed more 
transparently in the present formulation. I note that the 
2 J particles are coupled to each other, not only dynam- 
ically through the Hamilton function, but also kinemat- 
ically via the symplectic form. Eq. (|16p is the quantum 
mechanical counterpart of the Landau-Lifshitz equation 
for spin dynamics 20], which we would recover if we 
would restrict our description to (quasi-classical) CSs. 
For quantum spin systems, such as molecular magnets 
21j . the latter is clearly inadequate, and a fully quan- 
tum description as given in Eq. (fTB"|) is necessary. 

Finally, I briefly address the question of systems of 
interacting spins in magnetically ordered systems. The 



usual treatment, spin- wave theory, amounts to use a vari- 
ational wave-function given as a tensorial product of CS 
with site-dependent unit vectors n(R,-,t), and solve the 
linearized coupled Landau-Lifshitz equations. This ap- 
proach is clearly not adequate for systems with exotic 
(e.g., quadrupole or higher- multipole) ordering such as 
spin nematics 22[ , and the effective field theory proposed 
for spin nematics [23| applies only to spin-1 systems and 
cannot be extended to systems with spin J > 1. This 
clearly calls for a more general theory of spin nematics. I 
argue here that the most natural description of such sys- 
tems is the geometric description offered by the Majorana 
representation. It consists in setting up a quantum field 
theory based upon Majorana's constellations instead of 
the quasi-classical CSs; one thus obtains a path-integral 
theory for the 2 J coupled 0(3) fields u^. This theory will 
be developed in detail in a forthcoming paper. 

I gratefully acknowledge helpful discussions and/or 
correspondence with Robert Whitney, Efim Kats, and 
John Hannay. 



[6] 

m 

[8] 
[9] 
[10] 



[11] 
[12] 

[13] 
[14] 



[15] 



patrick.bruno@esrf.fr 

M.V. Berry, Proc. Roy. Soc. A 392, 45 (1984). 

Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 

(1987). 

A. Shapere and F. Wilczek (Eds.), Geometric phases in 
physics, World Scientific, Singapore (1989). 
A. Bohm et at, The Geometric Phase in Quantum Sys- 
tems, Springer- Verlag, Berlin (2003). 

D. Chruciriski and A. Jamiolkowski, Geometric Phases in 
Classical and Quantum Mechanics, Birkhauser, Boston 
(2004). 

P. Bruno, Phys. Rev. Lett. 93, 247202 (2004). 
J.H. Hannay, J. Phys. A 31, L53 (1998). 

E. Majorana, Nuovo Cimento 9, 43 (1932). 

F. Bloch and I.I. Rabi, Rev. Mod. Phys. 17, 237 (1945). 
J. Schwinger, US Atomic Energy Commission, Report 
NYO-3071 (1952); later published in Quantum Theory 
of Angular Momentum, ed. by L.C. Biendenharn and 
H. Van Dam, Academic Press, New York (1965). 

A. Perelomov, Generalized Coherent States and Their 
Applications, Springer- Verlag, Berlin (1986). 
J. Schwinger, unpublished (1945, 1959); reprinted in A 
quantum legacy: seminal papers of Julian Schwinger, ed. 
by K.A. Milton, World Scientific, Singapore (2000). 
J.M. Caillol, J. Physique Lett. 42, L245 (1981). Note that 
the expression of the 2D spherical Coulomb interaction 
was given there without proof; a proof may be found in 
Log-Gases and Random Matrices, P.J. Forrester, Prince- 
ton University Press (2010). 

Such an identification of a quantum state with the ther- 
modynamic partition function of a fictitious classical gas 
is familiar from the theories of the fractional quantum 
Hall effect (R.B. Laughlin, Phys. Rev. Lett. 50, 1395 
(1983)) and of the valence bond solids in quantum antifer- 
romagnets (D.P. Arovas, A. Auerbach and F.D.M. Hal- 
dane, Phys. Rev. Lett. 60, 531 (1988)). 
C.T. Lee, J. Phys. A 21, 3749 (1988). 



■5 



[16] R. Gilmore, J. Phys. A 9, L65 (1976). 
[17] J. Anandan and Y. Aharonov, Phys. Rev. Lett. 65, 1697 
(1990). 

[18] V.I. Arnold, Mathematical Methods of Classical Mechan- 
ics, Springer- Verlag, Berlin (1978). 

[19] P. Lebceuf, J. Phys. A 24, 4575 (1991). 

[20] L.D. Landau and E.M. Lifshitz, Phys. Z. Sowjetunion 
8, 153 (1935); reprinted in Collected Papers of Landau, 
edited by D. ter Haar (Pergamon Press, Oxford 1965), p. 
101. 



[21] D. Gatteschi, R. Sessoli, and J. Villain, Molecular Nano- 

magnets, Oxford University Press (2006). 
[22] H.H. Chen and P.M. Levy, Phys. Rev. Lett. 27, 1383 

(1971); A.F. Andreev and LA. Grishchuk, Sov. Phys. 

JETP 60, 267 (1984); P. Chandra and P. Coleman, Phys. 

Rev. Lett. 66, 100 (1991). 
[23] B.A. Ivanov and A.K. Kolezhuk, Phys. Rev. B 68, 052401 

(2003). 



